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Abstract 

The probability of evaporation induced by R + and R~ rotons at the surface 
of superfluid helium is calculated using time dependent density functional 
theory. We consider excitation energies and incident angles such that phonons 
do not take part in the scattering process. We predict sizable evaporation 
rates, which originate entirely from quantum effects. Results for the atomic 
reflectivity and for the probability of the roton change-mode reflection are 
also presented. 

PACS number: 67.40.-w, 67.40.Db 
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Quantum evaporation occurs in superfluid 4 He when a high-energy phonon or roton 
propagates to the surface where it annihilates and an atom is ejected in the free space (see, 
for example, Ref. This phenomenon is especially interesting because of the peculiar 

dispersion law exhibited by rotons. 

Despite the significant experimental and theoretical P-TH[ efforts made in the last 



years, the fundamental mechanisms underlaying the phenomenon of quantum evaporation 
are not yet understood. The experiments by Wyatt and co-workers |l]J^,§| have revealed 
that the main process is a one to one process (one excitation to one atom). This behavior is 
confirmed by the separate conservation of the energy and of the momentum parallel to the 
surface in the evaporation process. In contrast measurements of atom condensation |15j have 



instead pointed out the importance of non linear processes associated with the excitation of 
ripplons. 

The theoretical studies have not yet provided a clear and consistent picture of quantum 
evaporation. The reason is that it is very difficult to develop a reliable description of this 
phenomenon on a microscopic basis. In fact a good theory should be able to account for 
several effects: 

• (i) a correct description of the structure of the free surface, as well as of the elementary 
excitations of the system; 

• (ii) a quantum description of the scattering processes involving the elementary excita- 
tions at the surface; 

• (iii) the inclusion of inelastic channels (multi-phonons, multi-ripplons). 

A useful discussion concerning the role of quantum effects has been recently made in 
Refs. [0,0 where it has been pointed out that, due to the peculiar form of the maxon- roton 
dispersion exhibited by superfluid helium, there are severe constraints on the structure of the 
classical orbits associated with the elementary excitations when they cross the interface. In 
particular one finds that only phonons and rotons above the maxon energy (about 14 K) can 
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give rise to evaporation. Vice-versa, the theory of classical orbits predicts no evaporation 
from rotons with energy smaller than the maxon energy, because of the occurrence of a 
barrier at the interface. The experimental evidence for quantum evaporation induced 
by R + rotons, of all wave vectors, is consequently an important proof of the crucial role 
played by quantum effects. The quantum states associated with the above classical orbits 
(WKB states) have been also used to carry out a perturbative description of the scattering 
process ||14j| . However, the perturbative approach is not easily justifiable in this context. 
The above discussion has also an important implication on the reflection coefficient. In 
fact, according to the classical picture of Ref. ||13|| , atoms travelling at incidence angles 
sufficiently large, so that they lie in the shaded part of the energy diagram in Fig. p], are 
reflected with unit probability. (Note that in Fig. [I] and in the following, q x is the component 
of wave vector parallel to the surface, while z is the orthogonal coordinate). This behavior 
is contradicted by the experiments of Ref. ||15|| , which instead indicate full condensation 
also in that region. Actually experimental data show significant reflection only when the 
perpendicular component q z of the atom wave vector is close to zero ]TB| . 



The purpose of this Letter is to provide a first calculation of the evaporation rates using 
a many body approach accounting for both the requirements (i) and (ii) discussed above. 
We have chosen the values of q x and of the energy of elementary excitations in such a way 
that phonons do not take part in the scattering process; this corresponds to considering the 
shaded region in Fig. [l|. According to classical theory this region is characterized by full 
reflection and no evaporation. It is consequently an ideal region to test the importance of 
quantum effects. 

Our approach is based on the study of a slab of liquid 4 He confined in a box of size Lb ox , 
as shown in Fig. 0(a). The system is assumed to be translationally invariant in the x, y 
direction. The slab is chosen enough thick (50 ~ 70 A) to provide a quantitatively correct 
description of the behavior of the semi-infinite medium. The box size [Lbox > 100 A) has 
been chosen in order to allow a few oscillations of the atom wave function in the free space. 

We have calculated the eigenenergies and eigenfunctions of this system in the framework 



of density functional theory (see, for instance, Refs. [|I7,|I8| and references therein). When 
applied to a Bose system this theory describes the fluctuations of the density p and of the 
velocity potential </> according to the expansion 



= V>o(z) + f{z)e i{ -^ x -^ + g(z)e i(q * x+ujt) (1) 



where ipo{z) = y po(z) is fixed by the ground state density of the system and f(z) and 
g(z) are real wave functions to be determined, together with the frequency u, by solving 
self-consistently the equations of motion 

5 Jdtjdr (nty*, ip] - ipHh^il^J = (2) 

linearized with respect to / and g. The quantity E = Jdr 7i[ip* , ip] is the energy functional 
of the system (depending on p and 4>) which is assumed to be known. The same functional 
provides, through the variational procedure S(E — pN) = 0, the ground state profile po(z). 
In this work we use the density functional recently proposed in Ref. NTS]]. It provides a 
correct description of the equation of state of superfluid helium, as well as of the density 
profile at the surface. Furthermore it reproduces the dispersion law of phonons, maxons 
and rotons. The same theory gives a reliable description of surface modes (ripplons) both 
at small and high momenta fll7[ . 

Equations (|l|J^) have the typical form of the equations of the random phase approxima- 
tion (RPA). In particular they account for both particle-hole (f(z)) and hole-particle (g(z)) 
transitions which are coupled by the equations of motion (0). This coupling is of crucial 
importance in order to treat the correlation effects associated with the propagation of ele- 
mentary excitations in an interacting system. The equations of motion have also a structure 
formally identical to the one of the Bogoliubov equations for the dilute Bose gas and to the 
one of the Beliaev equations for Bose superfluids In the vacuum they coincide with 



the Schrodinger equation for the free atom wave function f(z), while g(z) vanishes. 

One should note that the equations of time dependent density functional (TDDF) theory 



correspond to fully quantum mechanical equations and consequently account for the interfer- 
ence and tunneling phenomena which are expected to play a crucial role in the evaporation 
process. Of course, due to linearization, they do not include inelastic processes associated 
with multi-phonons or multi-ripplons. These effects lie beyond the present formulation of 
the theory which is essentially a mean field theory. For the same reason in our theory 
phonons have an infinite life time and cannot decay into multi-phonons as is instead known 
to occur at energies below 10 K. Despite the absence of inelastic processes, we think that the 
solution of the evaporation problem within linearized TDDF theory is nevertheless highly 
instructive. 

The solution of the equations (0) can be determined with high precision working in the 
slab geometry discussed above. The solutions are real and either symmetric or antisymmetric 
with respect to the center of the slab. A typical solution is shown in Fig. 0(a), where we 
plot the function f(z), for an excitation at q x = 0.7 A -1 and %lo = 11 K. The figure shows 
the existence of atoms travelling outside the slab and of elementary excitations inside the 
slab. The corresponding function g(z), not shown, has also an oscillatory behavior inside the 
slab, while it vanishes outside consistently with the fact that the hole-particle components 
of the wave function ([I]) (associated with correlation effects) are absent in the free atom 
region. The wave length of the atom wave function is easily calculated starting from the 
energy conservation law Tiuj = [h 2 q 2 / (2m) — /i] where q 2 = q 2 + q 2 , while [i = —7.15K is the 
chemical potential of helium atoms. One finds X z = 2ir/q z = 16.4 A in agreement with the 
numerical results shown in Fig. 0(a). 

Due to the value of q x and u the solution shown in Fig. 0(a) cannot contain phonon 
components. This is best illustrated in Fig. 0(b) where we show the Fourier transform of 
the signal inside the slab. The signal reveals two distinct peaks, one corresponding to a R~~ 
roton with q z = 1.43 A -1 (q = 1.59 A -1 ), and the other to a R + roton with q z = 2.05 A -1 
(q = 2.17 A -1 ). In this case the R~ and R + rotons propagate at angles of 26° and 19° 
relative to the z-axis, respectively, while the atom outside the slab an angle 

of 61°. The additional oscillatory structure revealed by Fig. 0(b) originates from the finite 
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size of the box, smaller than L s z a &, used to calculate the Fourier transform. In fact one can 
well fit the calculated signal starting from a function of the form 

f{z) = f+ cos(g+z) + /_ cos{q~z) (3) 

inside the slab and f\z) = outside, as shown in Fig. 0(b). Actually this fitting procedure 
has been used in order to extract the values of /+ and /_ needed for the analysis of the 
evaporation rates. The same procedure has been used to analyze the function g(z). 

The results for / and g can be used to calculate the current giving the number of 
elementary excitations crossing the unit surface per unit time through the following equation: 

k = Vi(| fi\ 2 ~\ 9i T) (4) 

where {\fi\ 2 — \g%\ 2 ) is the density of elementary excitations and Vj is the group velocity of 
the i-th excitation (i = a, +, — ). The structure of the current (f|) emphasizes a remarkable 
feature of the RPA (or Bogoliubov) equations. Note that only for a free atom Eq. (Q) takes 
the familiar form j = \f\ 2 hc[/m. In a correlated system, like superfluid helium, the current 
behaves quite differently. For instance, in the long wave length phonon regime [q — ► 0), the 
group velocity coincides with the sound velocity and g — f ■ 

The current (||) is used to calculate the (real) amplitudes A4 = J\jf\sga(fi) of the signal 
relative to the various components (atom, i? ± rotons) in the scattering process taking place 
at the surface. The amplitude A a of the signal relative to the outcoming atom can be related 
to the ones of the incoming atom and R^ 1 rotons through the relation 

_ Aafr ntU = Saa A a e^ L « + iS^ a A_e-^ L 

+ iS +a A + e^ L (5) 

where Sij is the scattering matrix, 2L = L s i a b is the slab thickness and L a = (Lb ox — L s i a b)/2. 
Relation (|J) holds for symmetric states; a similar relation can be written for antysymmetric 
states. Notice that the phonon contribution to the scattering process is absent due the 
choice made for q x and for the energy. 
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The matrix elements of the scattering matrix entering Eq. ([5]) satisfy the relation SV,- = Sji 
which follows from unitarity and time reversal symmetry conditions. In terms of these matrix 
elements the evaporation probabilities P + and P~ (relative to R + and R~ rotons) and the 
reflection coefficient R take the form 

P + =| S +a | 2 ; P =| S- a | 2 ; R=\ S aa | 2 . (6) 

Furthermore, due to unitarity, one has P + + P~ + R — 1 (this is true only if one ignores 
inelastic channels, as in the present theory). 

In order to extract the physical coefficients P + , P~ and R it is necessary to obtain 
various solutions at the same values of q x and energy, involving different combinations of 
the atom and roton signals. This has been achieved by varying the thickness of the slab 
and the size of the box. Of particular importance, for our analysis, is the occurrence of 
the so called "resonance" states which are characterized by the absence of the atom signal 
(A a = 0) outside the slab, due to destructive interference. These states are useful because, 
as Eqs. (^) and (^) clearly show, they allow one to identify the ratio P + / P~ with |yl_/A + | 2 . 
Results are shown in Fie;. E¥ function of the roton energy for two values of q x (0.7 A 1 

and 0.8 A -1 ). The error bars are mainly due to the fact that the atom wave function is not 
exactly vanishing, even for the best resonant states resulting from the numerical solution 
of the equations of motion. This produces a statistical uncertainty in the values of the 
branching ratio. 

Clearly the determination of the reflection coefficient, as well as of the other elements 
of the scattering matrix, requires the analysis of non resonant states. A more detailed 
description of the procedure and a systematic discussion of the results will be presented in a 
longer paper. In Fig. |3](b) we show our results for the evaporation probabilities and for the 
reflection coefficient as a function of energy. The numerical uncertainty on these values is 
expected to be less than 10%. The analysis allows us to estimate also the probability P H 
for the roton change-mode reflection, R + «-> R~. We obtain P + ~ ~ 0.3 and 0.2 at %u = 10.8 
K and 11.3 K, respectively. 
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The main conclusions emerging from our results are: 

• Quantum effects give rise to sizable evaporation rates of rotons in the region of energy 
and angles where evaporation is not allowed classically. 

• R~ rotons turn out to be less active in the evaporation process than R + rotons. Evap- 
oration from R + rotons becomes dominant when the energy increases. 

• The probability for the roton change-mode reflection is sizable in the energy interval 
considered and decreases with energy. 

• The atom reflection coefficient is smaller than 10% for energy greater than about 11 
K, and decreases for higher energies. 

The value of the reflection coefficient below 11 K is still too large with respect to the 
experimental data, but nevertheless its sizable decrease from the classical value R — 1 
reveals the very important role played by quantum effects. The remaining discrepancy with 
experiments is likely associated with inelastic processes not accounted for in the present 
calculation. 

The above results concern the region of large q x and large incident angles, where phonons 
do not take part in the process. For normal impact and energy smaller than the roton 
minimum (but larger than 7.15 K), one has the opposite regime, where only phonons take 
part in the scattering process. The results of our calculations in this case give very small 
values for the atom reflection coefficient, in agreement with experiments. 

We are indebted to C. Carraro and A.F.G. Wyatt for many fruitful discussions. This 
work was partially supported by the US Department of Energy Office of Basic Sciences under 
contract no. W-31-109-ENG-38. 
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FIGURES 

FIG. 1. Spectrum of elementary excitations in superfluid 4 He. Solid line: phonon-maxon-roton 
dispersion; dashed line: threshold for atom evaporation; dot-dashed line: dispersion of the free 
surface mode; shaded area: region of roton-atom processes. 

FIG. 2. An example of solution f(z) for q x = 0.7A -1 , hu; = UK, L s[ab = 62A and L box = 140A. 
The dashed line in the upper part (a) is the density profile of the slab, in arbitrary units. In the 
lower part (b) the Fourier transform of f(z) inside the slab is also shown. The best fit with formula 
@ (dashed line) is practically indistinguishable from the numerical solution (solid line). 

FIG. 3. Ratio of the P + and P~ evaporation probabilities (top) and absolute values of the 
evaporation and reflection probabilities (bottom) as a function of energy. Triangles, circles and 
squares correspond to P + , P~ , and R, respectively. All values below 11. 5K are calculated at fixed 
parallel wave vector q x = 0.7 A -1 , the others at q x = 0.8 A -1 . 
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